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Multi-particle production in QCD is dominated by higher twist contributions. The opera- 
tor product expansion is not very effective here because the number of relevant operators 
grow rapidly with increasing twist. The Color Glass Condensate (CGC) provides a frame- 
work in QCD to systematically discuss "classical" (multiple scattering) and "quantum" evolu- 
- - ■ tion (shadowing) effects in multi-particle production. The apparently insuperable problem of 

' nucleus-nucleus scattering in QCD simplifies greatly in the CGC. A few examples are discussed 

, with emphasis on open problems. 

(N 

Q ' 1 Introduction: the twist expansion and small x physics 

In the Bjorken limit of QCD, oo, s — > oo, xbj ~ Q'^/s = fixed, we have a powerful 

framework to compute a large number of processes to high accuracy. Underlying this machinery 
is the Operator Product Expansion (OPE), where cross-sections are identified as a convolution 
i of short distance "coefficient functions" which are process dependent and long distance parton 

I distribution functions which are universal. The evolution of the parton distribution functions 

with X and is described by splitting functions, which determine the probability of "parent" 
partons to split into a pair of "daughter" partons. Both coefficient functions and splitting 
functions for DIS inclusive cross-sections are now available to Next-Next-Leading-Order (NNLO) 
accuracy . 

The leading contributions in the OPE come from operators which have minimal twist, where 
twist is defined as the dimension minus the spin of the operators. Higher twist operators are 
suppressed by powers of (m-^/Q-^)-^"-^™'", where m is a hadronic mass scale, T denotes the twist 
of the operator and Tmin the minimal twist (2 for the structure functions Fi and F2). These 
. higher twist operators can therefore be ignored in the Bjorken limit, albeit their contribution 

provides a systematic error to the application of the leading twist formalism at finite Q^. 

However, the bulk of multi-particle scattering in QCD is dominated by soft and semi-hard 
physics. In the language of the OPE, the latter are higher twist effects. These are of two 
5^ , sorts. The first are "kinematic" higher twist contributions, which arise from trace contributions 

that are often ign ored in the OPE, where the leading contributions are from symmetric and 
traceless operators^^EI, These kinematic high twist contributions are of order x^-^vn? /Q"^ and are 
of decreasing importance at small xbj- The other "dynamical" higher twist contributions are 
from the hadronic matrix elements of the higher twist operators themselves. The relevant twist 
four matrix elements for leptoproduction were discussed in great detail by Ellis, Furmanski and 
Petronzio'^ and by Jaffe and Soldate^^, and expressed in terms of multi-parton distributions 
by Jaffe ^. Discussions specific to twist four contributions at small x are contained in Ref. ^ 
There are many more contributions at twists greater than four-these have not been quantified. 

To understand why the twist expansion is likely not a useful expansion at small x, we need to 
consider the properties of the theory in the Regge limitrxBj ^ 0, s — > cxd, = fixed. The BFKL 
renormalization group equation'^ describes the leading as ln(l/x) behavior of gluon distributions 
in this limit. The solutions of the BFKL equation predict that gluon distributions grow very 
rapidly with decreasing x. In the Regge asymptotics, since the transverse size of the partons is 
fixed, this growth of distributions will lead to the overlapping of partons in the transverse plane 
of the hadron. In this regime, contributions that were power suppressed in the BFKL scheme 
become important. These are recombination and screening effects which slo w dow n the growth 
of gluon distributions leading ultimately to a saturation of these distributionJ^^^I. Such effects 
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must appear at small x because the occupation number" of partons in QCD can be at most of 
order l/as- 

Qualitatively, the competition between Bremsstrahlung and recombination/screening ef- 
fects become of the same order when ^ ^^j^iQ'i ^ ~ as{Q'^) ' "^^^^ relation is solved self- 
consistently when Q = Qs{x). The scale Qsix) is termed the saturation scale and it grows as 
one goes to smaller values of x. What does this have to do with higher twists? As we will 
discuss further in the next section, when Qsix)"^ > Q^, all higher twists are equally important. 
The OPE therefore is not a good expansion in this small x kinematic region, where the typical 
momentum of partons is of order Qg ^. 

There is however a glimmer of hope in this seemingly hopeless situation. This is because 
Qlix) » Aqqjj, which suggests that weak coupling techniques in QCD are applicable in the 
Regge limit. In the next section, we will discuss a weak coupling effective field theory approach 
which provides a more efficient organizing principle than the OPE at high parton densities. 



2 The Color Glass Condensate 

The physics of high parton densities can be formulated as a classical effective theory 
because there is a Born-Oppenheimer separation between large x and small x modes'^ which 
are respectively the slow and fast modes in the effective theory. Large x partons are static 
sources of color charge for the dynamical wee (small x) parton fields. The generating functional 
of wee partons has the form 

where the wee parton action has the form 

S[A,p] = ^ J d^xFl, + j^ j d^x^dx-5{x-)Ti{p{x^)U,oo,oo[^-]) • (2) 

In Eq. n /3 is a two dimensional classical color charge density and W[p\ is a weight functional 
of sources (which sit at momenta fc+ > A+: note, x = ^"""/-fhtidron)- "^^^ sources are coupled to 
the dynamical wee gluon fields (which in turn sit at k'^ < A+) via the gauge invariant term'^ 
which is the second term on the RHS of Eq. |2j Here C/-oo,oo denotes a path ordered exponential 
of the gauge field A~ in the x+ direction. The first term in Eq. [21 is the QCD field strength 
tensor squared-thus the wee gluons are treated in full generality in this effective theory, which is 
formulated in the light cone gauge = 0. The source j is an external source-derivatives taken 
with respect to this source (with the source then put to zero) generate correlation functions in 
the usual fashion. 

The argument for why the sources are classical is subtle and follows from a coarse graining of 
the ef fective action. The weight functional for a large nucleus is a Gaussian in the source den sity 
| 14 | 1^ |^ with a small correction for SU(iVc) coming from the Nc — 2 higher Casimir operators'^. 
The variance of the Gaussian, the color charge squared per unit area proportional to is 
a large scale-and is the only scale in the effective action*^. Thus for Aq(jj-), a5(^^) << 1, 

"This corresponds to the number of partons per unit transverse area, per unit transverse momentum, in light 
cone gauge. This condition has its gauge invariant counterpart in the requirement that the field strength squared 
not exceed l/as. 

Another more severe reason why th e O PE breaks down-even at leading twist at small x-has to do with the 
infrared diffusion in the BFKL equation^! Ironically, this diffusion is cured by higher twist satu ration effects. 
■^For an alternative gauge invariant form, which also recovers the BFKL equation, see Ref.'^. 

is simply related in the classical theory to the saturation scale via the relation Ql = 
asNcA ln(Qs/A|cD) 



and one can compute the properties of the theory in Eq-^in weak coupHng. For an SU(3) Yang- 
Mills theory, there is an additional contribution proportional to the cubic Casimir operator. It 
is parametrically suppressed by A ^/^- However, this term generates Odderon excitations in the 
CGC already at the classical level^^. 

The saddle point of the action in Eq. Ogives the classical distribution of gluons in the nucleus. 

The Yang-Mills equations can be solved analytically to obtain the classical field of the nucleus 

1/1 1 f-1 on 

as a function of p: ^ci.(/o) ' " ■ One can determine, for Gaussian sources, the occupation 
number cp = dN /-kE? / dk\dy of wee partons in the classical field of the nucleus. One finds for 
k±_ » the Weizsacker- Williams spectrum ~ Ql/k\\ for k± < Qs, one obtains a complete 
resummation to all orders in k±, which gives ~ ^ ^T^{Qs/k±). (The behavior at low k± can, 
more accurately, be represented as — r(0, z) where F is the incomplete Gamma function and 
z = A^i/Q^O). 

Small fluctuations about the effective action in Eq. |21 give large corrections of order 
a5ln(l/a;). The Gaussian weight functional is thus fragile under quantum evolution of the 
sour ces. A Wilsonian renormalization group (RG) approach systematically treats these correc- 
tionsl^. In particular, the change of the weight functional W[p] with x is described by the 
JIMWLK- non-linear RG equations'^. These equations form an infinite hierarchy of ordinary 
differential equations for the gluon correlators < A1A2 ■ ■ ■ An >y, where Y = ln(l/x) is the ra- 
pidity. For the gluon density, which is proportional to a two-point function < a°'{x±)a''{y±) >, 
one recovers the BFKL equation in the limit of low parton densities. Further dev elopments 
beyond the JIMWLK equation have been summarized at this conference by lancu^^. 

In the limit of large Nc and large A {ogA^/^ >> 1), the JIMWLK hierarchy closes for the two 
point correlator of Wilson lines because the expectation value of the product of traces of Wilson 
lines factorizes into the product of the expectation values of the traces: {Tt(VxV^)Tt{VzV^)) — > 
{Tt{V^V^)) (TriVzV^)), where K = VexpiJ dz' a" {z' , x ^)T'') . Here V denotes path ordering 
in x~ and is an adjoint SU(3) generator. The cross-section for a qq pair scattering off a target 
can be expressed in terms of these 2-point dipole operators as (JqqN^x, r±) = 2 J d^h My{x, r±, b), 
where My, the imaginary part of the forward scattering amplitude, is defined to be A/y = 
1 — ^ < Tr{VxV^) >Y- The size of the dipole, f± = x± — y± and b = {x± + y±)/2. The 
JIMWLK equation for the two point Wilson correlator is identical in the large A, large Nc 
mean field limit to a n eq uation derived independently by Balitsky and Kovchegov-the Balitsky- 
Kovchegov equation I^HI^ which has the operator form = ag/CsFKL ® {My — My}- Here 
/Cbfkl is the well known BFKL kernel. When AA << 1, the quadratic term is negligble and 
one has BFKL growth of the number of dipoles; when M is close to unity, the growth saturates. 
The B-K equation is the simplest equation including both the Bremsstrahlung responsible for 
the rapid growth of amplitudes at small x as well as the repulsive many body effects that lead 
to a saturation of this growth. 

We now return to our discussion of higher twists in the previous section. In lepto-production, 
the structure function F2 at small x is proportional to \^p-y*^gg\'^ ®crqqNi where I'i/'^*^^^^ is the 
probability for a virtual photon to split into a qq pair and aqqN is the qq — N cross-section 
discussed previously. Since the latter is proportional to a product of Wilson lines, F2 gets 
contributions from A^-point gluon distributions. In the classical McLerran-Venugopalan (MV) 
model of Gaussian color sources, these can be expressed explicitly!^ as an expansion in Q^/Q^- 
thus, for > Q^, all higher twists contribute equally. The OPE would not be very useful in 
this region-however, in the CGC framework, higher twist effects are included both at the tree 
level in the MV model, and in the small x quantum evolution of the BK and JIMWLK RG 
equations. It is interesting to ask whether the effective theory of the CGC at sufficiently small 
X and large can be matched on to the full theory beyond the leading twist level-this has not 
been done thus far but is feasible in principle. 



3 Hadronic Scattering in the CGC 



At small X, both the collinear factorization and k± factorization limits of pQCD can be un- 
derstood in a systematic way in the framework of the CGC. Rather than a convolution of 
probabilities, one has instead the collision of classical gauge fields. The expectation value of an 
operator O can be computed as 

<0>Y=J [dpi] [dp2] W,, [pi] W,, [p2] 0{pi,p2) , (3) 

where Y = \n{l/xF) and xp = xi — X2- All operators at small x can be computed in the 
background classical field of the nucleus at small x. Quantum information, to leading logarithms 
in X, is contained in the source functionals Wxj^(x2)[PiiP2)]- The operator O can be expressed in 
terms of gauge fields A^^[pi, p2]{x). 

Inclusive gluon production in the CGC is computed by solving the Yang-Mills equations 
[D^,F'^'^]"' = J^'°-, where J'^ = pi 6{x~)6'^^ + p2d{x~^)6'^~ , with initial conditions given by the 
Yang-Mills fields of the two nuclei before the collision. These are ob tain ed self-consistently by 
matching the solutions of the Yang-Mills equations on the light cone^^. The initial conditions 
are determined by requiring that singular terms in the matching vanish. Since the Yang-Mills 
fields in the nuclei before the collision are known, the classical problem is in principle com- 
pletely solvable. Quantum corrections not enhanced by powers of Q5ln(l/x) can be included 
systematically. The terms so enhanced are absorbed into the weight functionals TV[/3i,2]- Thus 
all "classical" multiple scattering effects are obtained by solving the Yang-Mills equations, while 
the small x quantum evolution effects (which gives rise to shadowing) are contained in the weight 
functionals which obey the JIMWLK/BK equations. 

Hadronic scattering in the CGC can therefore be studied through a systematic power count- 
ing in the density of sources in powers of pi,2/^i^-i2- The power counting is applicable either 
to a proton at small x, or to a nucleus (whose parton density at high energies is enhanced by 
A^/'^) at large transverse momenta. The relevant quantity here is Qs, which, as one may recall, 
is enhanced both for large A and small x. So as long as k± >> Qg » Aqcd, the proton or 
nucleus is considered dilute. One can begin to study the applicability of both collinear and k± 
factorization at small x in this approach. 

To lowest order in ppi/k\ and Pp2/k']_, one can compute inclusive gluon production ana- 
lytically At large transverse momenta, Qs « k±, the scattering can be expressed in a 
/c_L-factorized form. The inclusive cross-section is expressed as the product of two unintegrated 
{k± dependent) distributions times the matrix element for the scattering. k± factorization is a 
good assumption at large momenta for quark pair-production. This was worked out in the CGC 
appr oach by Frangois Gelis and myself In this limit, our result agrees exactly with that of 
Refs.EHI. 

In the semi-dense/pA case, one solves the Yang-Mills equations to determine the gluon field 
produced-to lowest order in the proton source density and to all orders in the nuclear source 
dens ity. Th e inclusive gluon production cross-sect ion , in this fra mework, was first computed in 
j^gfg lsii]l^ and shown to be k± factorizable in Ref.^^. In Ref.^^, the gluon field produced in pA 
collisions was computed explicitly in Lorentz gauge d^A^ = 0. The fact that the distributions 
are k± factorizable is remarkable because the "unintegrated" gluon distribution of the nucleus is 
not the usual leading twist unintegrated distribution, but includes all higher twists. Its evolution 
with energy is given by the JIMWLK/BK equations. 

If we wish to study multiple scattering effects alone, these can be studied in the MV model- 
which provides the initial conditions for quantum evolution. Thus at larger x's multiple scat- 
tering effects dominate while quantum effects turn on as one goes to smaller x due to the RG 
evolution of the weight functionals. The well known "Cronin" effect is obtained in our formalism 



and can be simply understood in terms of the multiple scattering of a parton from the projectile 
with those in the target. The remarkable energy dependence of the Cronin effect seen by the 
RHIC experiments may be due to quantum evolution effects and has been discussed elsewhere^^l. 

Quark production in p/D-A collisions can be computed with the gauge field in Lorentz 
gauge^^. Unlike gluon production, neither quark pair-production nor single quark production 
is strictly k± factorizable. The pair production cross-section however can still be written in k± 
factorized form as a product of the unintegrated gluon distribution in the proton times a sum of 
terms with three unintegrated distributions, 4'g,g^ ^qq,g ^ggm- These are respectively pro- 
portional to 2-point, 3-point and 4-point correlators of the Wilson lines we discussed previously. 
For instance, the distribution 4>qq,g is the product of fundamental Wilson lines coupled to a qq 
pair in the amplitude and adjoint Wilson lines coupled to a gluon in the complex conjugate 
amplitude. For large transverse momenta or large mass pairs, the 3-point and 4-point distribu- 
tions collapse to the unintegrated gluon distribution, and we recover the previously discussed 
fej^-factorized result for pair production in the dilute/pp-limit. Single quark distributions are 
straightforwardly obtained and depend only on the 2-point quark and gluon correlators and 
the 3-point correlators-which are "all twist" operators as previously. For Gaussian sources, as 
in the MV-model, these 2-, 3- and 4-point functions can be computed exactly as discussed in 
Ref . The renormalization group evolution of these distributions lead to shadowing of the 
distributions. Understanding their evolution with energy may provide important information 
about the structure of multi-parton correlations in high energy QCD. 

Our results for gluon and quark production in p/D-A c ollis ions (for a review, see Ref.l^^, 
coupled with the previous results for inclusive and diffractivel^ distributions in DIS suggest an 
important new paradigm. At small x in DIS and hadron colliders, quark and gluon structure 
functions, which are the right observables in a leading twist formalism, are no longer the right 
observables to capture the relevant physics. Instead they should be replaced by dipole and 
multipole correlators of Wilson lines that seem ubiquitous in all high energy processes and are 
similarly gauge invariant*^ and process independent. For a similar conclusion on the importance 
of unintegrated distributions from a different perspective, see Ref."^. To determine whether 
these distributions are robust, next-to-leading order computations must be performed, which is 
a formidable, but by no means impossible, task. 

In nucleus-nucleus collisions, pi,2/^j_ ~ 1- There is no small expansion parameter and one 
has thus far not been able to compute particle production analytically in the C GC. Uii like gluon 
production in the pp and pA cases, A;_|_-factorization breaks down in the AA- 

cas JM. A si gnif- 

icant consequence is that one cannot factor the quantum evolution of the initial wavefunctions 
into unintegrated gluon distributions unlike the pA case. Nevertheless, there is a systematic 
way to include small x effects in the AA case. The problem of nuclear collisions is well defined 
in weak coupling and can be solved numerically"^^. The numerical simulations thus far assume 
Gaussian distributions of the color sources as in the MV model. This is reasonable for central 
Gold-Gold collisions at RHIC where the typical x is of order 10^^. At the LHC, the typical x at 
central rapidities is an order of magnitude lower. At these x values, quantum evolution effects 
are important and one should use one solutions of JIMWLK/BK RG equations to determine 
the distribution of sources. 

We will restrict ourselves to discussing numerical solutions for Gaussian color sources. The 
saturation scale Qs (which is an input in the numerical solutions in this approximation) and the 
nuclear radius R are the only parameters in the problem. The energy and number respectively 
of gluons released in a heavy ion collision of identical nuclei can therefore be simply expressed 
^ ^FH^ ^ ^^'^ ^ ~ ^ Qli "^here (up to 10% statistical uncertainity) we compute 

"^This relies on requiring that one closes the ends of the Wilson lines at ±cx) in the nucleus to form a closed 
loop, and further that the fields go to zero at ±oo. While physically attractive, this construction is arbitrary. I 
thank F. Gelis and K. Kajantie for a discussion on this point. 



numerically 1^ ce = 0.25 and cat = 0.3. Here r] is the space-time rapidity. The number 
distributions of gluons can also be computed in this approach. Remarkably, one finds that 
a) the number distribution is infrared finite, and b) the distribution is well fit by a massive 
Bose-Einstein distribution for k±/Qs < 1.5 GeV with a "temperature" of ~ 0A7Qs and by the 
perturbative distribution Q^/k^ for k±/Qs > 1.5. We will not discuss comparisons of CGC 
predictions to RHIC data here but refer to Ref. . 

The transition to the QGP from the CGC remains as an outstanding theoretical problem. 
Due to the rapid expansion of the system, the occupation number of modes falls well below one 
on time scales of order 1/Qs- From these times onwards, the classical approach breaks down 
for all but the softest modes-well before thermalization. On the other hand, for elliptic flow 
from hydrodynamics to be significant, the conventional wisdom is that thermalization should 
set in early. A necessary condition is that momentum distributions should be isotropic. The 
CGC initial conditions are very anisotropic with < p± >~ Qs and < pz >~ 0. How does this 
isotropization take place? All estimates of fina l state re-scattering of p artons formed from the 
melting CGC, both from 2^2 processes^^ and 2 — > 3 processes!^ suggest thermalization 
takes longer than what the RHIC collisions seem to suggest. In the "bottom up" scenario^^, 
Tthermai ' ^ A;, which at RHIC energies gives rthermai ~ 2 - 3 fm. 

Recently, it has been suggested that collective instabilities, analogous to the well known 
Weibel instabilities in plasma physics, can speed up themalization ■ . For a nice recent review, 
see Ref.l-^. Starting from very anisotropic (CGC-like) initial conditions, these instabilities drive 
the system to isotropy. In some estimates, these may be very short time scales of order 1/Qs- 
What is the relation of this language of instabilities and that of our classical field simulations? 
Paul Romatschke and I^^ recently looked at the effect of small violations of boost invariance 
on the dynamics of gauge fields produced with CGC initial conditions. We found that there is 
a Weibel instability and that the maximally unstable modes grow as exp(-y/r) (as opposed to 
exp(r) in a static box). This behavior was predicted in Ref.^^. Increasing the initial amplitude 
of the fields that violate boost invariance, we find that the growth of the instability terminates 
when non-Abelian effects become large. These studies don't address the possible isotropization 
of the system by instabilities because the initial amplitudes examined are much smaller than 
physically plausible. This work is in progress. 

An equally interesting problem is that of chemical equilibration. At high energies, the initial 
state in a heavy ion collision is dominated by gluons. Are quarks produced in sufficient numbers 
for the system to reach chemical equilibrium (where the ratio of gluons to quarks is expected 
to be 32/21Nf)? Recently, the Dirac equation for quarks was solved by Cells, Kajantie and 
Lappil^ in the background field computed in Ref. For as values comparable to those in 
the pre-equilibrium phase of RHIC collisions, they indeed find that sufficient numbers of quarks 
pairs are produced for chemical equilibrium to be plausible. As for the studies of collective 
instabilities, more intensive numerical studies are required for conclusive results. 
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